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Abstract. It is known that the imaginary parts of the frequencies of the quasi 
normal modes of the Schwarzschild black hole are equally spaced, with the level 
spacing dependent only on the surface gravity. We generalize this result to a wider 
class of spacetimes and provide a simple derivation of the imaginary parts of the 
frequencies. The analysis shows that the result is closely linked to the thermal 
nature of horizons and arises from the exponential redshift of the wave modes 
close to the horizon. 



In the study of small perturbations around a black hole spacetimes, quasi-normal 
modes (QNM, hereafter) play a crucial role and have been investigated very extensively 
(for a review, see It is known that the frequencies of QNM for the Schwarzschild 
black hole have the structure 



where k is the surface gravity of the black hole. This result was originally obtained 
numerically [3] and, surprisingly enough, was given an analytic proof only recently 
01 Ej. Most of the theoretical studies were either model dependent or concentrated on 
the real part of k n because of its relationship to loop quantum gravity and the area 
spectrum of black holes (see, for an incomplete sample of references, [HI El El EJ| ) • It 
is, however, equally intriguing why the imaginary part of k n has a, simple, equally 
spaced structure with the surface gravity of the horizon determining the level spacing. 
We shall try to provide a simple derivation of this feature in a context more general 
than that of a Schwarzschild black hole in this paper. [While this work was nearing 
completion, a paper |l(Jj appeared in the arXive with substantially same results. We 
shall provide brief comparisons wherever relevant.] 

For our analysis, we shall choose a general class of spherically symmetric metrics 
of the form 



It was shown in that spacetimes described by the above class of metrics have 
a fairly straight forward thermodynamic interpretation and — in fact — Einstein's 
equations can be expressed in the form of a thermodynamic relation TdS = dE — PdV 
for such spacetimes. Since quasi normal modes are conjectured to have possible, 
though as yet unclear, relationship with horizon thermodynamics, it will be interesting 
to explore the QNM for these metrics. We shall show, by a fairly straight forward and 
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simple analysis that the imaginary part of the QNM frequencies for these spacetimes 
is given by Im k n = inn for sufficiently large n. Our analysis shows that this result has 
the same mathematical origin as the standard Hawking evaporation, in the sense that 
both arise due to similar mathematical operations involving infinite redshift surface 
and tortoise coordinate. What is more, we shall only use the behaviour of f(r) on one 
side of the horizon; for example, in the case of Schwarzschild black hole we will obtain 
the results sticking to the r > 2M region, unlike for e.g., the results in ref. |Hj which 
depends on the form of the metric near r = 0. 

We are interested in the case in which the spacetime has one (or more) horizon(s) 
at which f(r) = and fir) is a finite constant. (The second condition rules out 
extremal black holes and similar spacetimes.) The surface gravity is then given by 
k = (1/2)| y' (a.) | where a is the radius of the horizon. For simplicity, we shall consider 
the single horizon situation for most of the analysis and will comment on generalization 
for multi horizon scenarios (which turns out to be far from trivial) in the end. It was 
shown in that the Einstein's equations allow us to write /(r) in the form 

f(r) = 1 - - - - /' p{r)r 2 dr (3) 



where p(r) is the energy density and we have chosen the normalization such that / = 
at the r — a surface which acts as a horizon. This provides an explicit context for our 
discussion. 

Consider a massless scalar field <j> satisfying the wave equation Dip = in this 
spacetime. (Generalization to higher spins are straight forward). We look for solutions 
to the wave equation in the form 

= IjP(r)r <m (0,0)exp(tfci), (4) 
r 

Straight forward algebra now leads to a "Schrodinger equation" for F given by: 
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F = k 2 F (5) 



where 



and r* is the tortoise coordinate defined by 

dr , . 

In equation (jSJ) we treat r* as the independent variable and all other quantities arc 
expressed in terms of r* through the function r»(r). Near the horizon r = a, we can 
expand f(r) in a Taylor series and use J7J) to obtain an explicit relationship between 
r and r„: 

1 lnf--l)+if--l N )+(!?[(r-a) 2 ] (8) 



2k \a ) \a 
where the length scale L depends on the source distribution p(r). [At a fundamental 
level, so does the surface gravity k; however, we shall see that our results are 
independent of L but depends on k.] With this definition, r > a is mapped to 
> — oo, and we will not require the region (r < a) beyond the horizon. 
We shall attempt to determine the frequencies of QNM by treating equation (J5J 
as a one-dimensional Schrodinger equation. In that case, the frequencies of QNM 
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can be identified with the poles of the scattering amplitude S(k) in the momentum 
space. The scattering amplitude, in turn, requires a complete solution of equation (J5J 
which, of course, is not available. We shall therefore invoke two approximations to 
obtain the scattering amplitude. First, we shall invoke the Born approximation for 
S(k), which should be valid for large momentum transfer and hence for n 3> 1 limit 
of QNM. Second, we shall evaluate the integral in Born approximation by taking the 
leading order contribution near the horizon. We shall first obtain the result and then 
comment on the validity of the approximations. 

The scattering amplitude in the Born approximation is given by the Fourier 
transform of the potential V(x) with respect to the momentum transfer q = kf — kj. 
In one dimension, k^ and kf should be parallel or anti-parallel; further we can take 
their magnitudes to be the same for scattering in a fixed potential. Then non-trivial 
momentum transfer occurs only for kf = — kj so that q = 2ki in magnitude. As we 
shall see, the factor 2 in this relation is crucial to get the correct result and it arises 
due to the constraint k, • + kf = 0. (Our original problem was three-dimensional and 
we are not working out the three-dimensional scattering amplitude in, say, s-wave 
limit. Rather, we first map the problem to an one-dimensional Schrodinger equation 
and study the scattering amplitude in one dimension). The scattering amplitude can 
now be expressed as 

/oo 
dr,V(r,)e 2lkr * (9) 
-oo 

where we have omitted irrelevant constant factors. Changing the integration variable 
to r using J7J and the form of the potential in J§J, we get 
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where e = [(r/a) — 1]. This integral, unfortunately, cannot be expressed in closed 
form (except in a few special cases; see equation (|15l) below). However, it is obvious 
that the dominant contribution, which is of interest to us, arises close to the horizon 
near e « 0. (This becomes clearer if we rotate the contour of integration to purely 
imaginary values of e by e — > it so that the integrand takes the form Q(x)x la e~ x 
with Q being regular near x = 0). Since /' is regular at e = 0, the leading order 
contribution near e = has the structure 

S(k) oc J™ ^e ik / K e i2ku oc T (l + i^j (11) 

where we have ignored a finite proportionality constant which is irrelevant for 
determining the poles. The poles of the scattering amplitude occur when the argument 
of the gamma function take negative integral values. This immediately gives the 
imaginary part of the QNM frequencies to be 

k n = inn (for n ^> 1) (12) 

This completes our derivation and we shall now comment on several features of the 
analysis. 

To begin with, let us note that there is a simple connection between the 
above result and the standard derivation of Bogoliubov coefficients in the case of 
a horizon which causes infinite redshift. A wave mode, the frequency of which 
is being exponentially redshifted, will be described by the amplitude of the form 
<j>(t) = exp(iqe~ Kt ) where q is a constant. Such a form arises in several contexts like 
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the black hole spacetime, Rindler frame and in a universe which is asymptotically De 
Sitter. An observer using the time coordinate t will decompose this mode into positive 
and negative frequencies with respect to t by the relation 

^ !{v)e ~ lvt (13) 

Elementary calculation shows that f{v) has the factors r(— iv/n) and T(l — iv/k) 
for positive and negative values of v; hence f(v) has poles exactly as in 1)12[). The 
integrals involved in this problem are identical to those we encountered above to 
obtain the scattering amplitude. The gamma function arises in both the cases due to 
the logarithmic behaviour of the tortoise coordinate or, equivalently, the exponential 
behaviour of r as a function of r* near the horizon. Thus, the structure of imaginary 
part of the QNM frequencies is closely linked to the standard analysis performed for 
deriving the thermal nature of horizons. 

In the case of Schwarzschild (or De Sitter spacetime) , the integral in (|10fl can be 
evaluated in closed form in terms of confluent hyper-geometric function \F% using the 
result: 

poo a — x 

I= dx 7T \ r =fr~"r(l + Q-n)iiq(n, -a + n,^ 1 ) 
Jo (1 + " x ) n 

+ W~Y ^ifil+a^ + a-^r 1 14 

r(n) 

In the case of Schwarzschild spacetime, / = 1 — ifl/r) and the scattering amplitude 
(omitting unimportant constant prefactors) becomes with 6 = k/n, 

, . s. 2-i8 

5(A) = f^J \T(1 + i6)\ 2 \ i l(l + 1)^(1 + id, i9,-2ikL) 

+ -(1 - i6) iFi(l + i6, id - 1, -2ikL) 
+ r(-2 + iO) [iFi(3, 3 - i0, -2ikL) 

- l ^r^-(2i + e) 1 F 1 (2,2-ie,-2ikL)] (15) 
ZkL J 

This exact expression confirms the structure of the poles identified earlier; the Gamma 

function has poles when its argument takes negative integral values and \F\(x, y,z) 

has poles when y takes negative integral values or zero. All these poles are summarized 

by equation l|12|) for large n. It is clear that no new poles appear at least in the case 

of Schwarzschild geometry when the exact analysis is performed. 

The integral in l|10(l can also be evaluated by retaining the full Taylor series 

expansion in e as is done in ref . |10j . Since each term of the power series is a gamma 

function, our key result is not changed. However, it is not clear whether such an 

expansion, retaining all the powers of e, is consistent with the Born approximation. 

It seems more natural to relate the result to the behaviour near the horizon which is 

what we have done. This is particularly true if the spacetime has a second horizon at 

r = b, say, in which case it is probably incorrect to extend the range of integration to 

all values of r. There is also the issue of new poles appearing when an infinite series 

is summed up which are not present in the sum containing finite number of terms 

(like the pole at x = 1 which arises when the series (1 + x + x 2 + • • •) is summed up 

retaining all the terms). This does not happen in the case of Schwarzschild geometry, 
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as is evident from the exact result in equation (|15f) and is unlikely to happen in a 
more general situation. We therefore can conclude that the poles of the scattering 
amplitude in the Born approximation occur along the imaginary axis. 

While the analysis appears to be strikingly simple, there are some subtleties which 
are not fully resolved. Since the result arises from the nature of the potential near 
the horizon, it seems reasonable to use an approximate form of the potential valid 
near the horizon in the equation J^} and solve for the scattering amplitude. It is fairly 
straightforward to do this (see, e.g., ref.pj])- To the lowest order, the horizon geometry 
can be approximated by that of a Rindler spacetime and the solutions to equation (JSJ 
can be expressed in terms of the Macdonald function K{ V {z). The scattering amplitude 
has unit modulus but exhibits divergent phase on the poles determined by equation 
Ijl2jl . At this order, only the surface gravity of the horizon plays any role and the 
analysis is equivalent to what we have done. Mathematically, this corresponds to 
retaining the lowest order contribution in the term in the square brackets in equation 
<|1U|I . The next order approximation to the potential near the horizon is provided by 
retaining one higher order in e in evaluating the integral in (|10() which will lead to 
the sum of two Gamma functions in our approach, without changing the structure of 
the poles. However, a different picture emerges if we solve the Schrodinger equation 
exactly in this approximate potential. When the potential is approximated by retaining 
one order higher than the Rindler case, then a more complicated form of scattering 
amplitude results (see equation 11 of ref.^21)- The imaginary part of the QNM 
frequencies obtained by this method has a level spacing which is 2k rather than the 
standard result k (see equation 12 of ref.^Hj; this is possibly related to the difference 
between q — 2ki and ki in the case of Born approximation mentioned earlier). It is 
not clear whether approximating the potential, to a finite order of accuracy near the 
horizon and then solving the Schrodinger equation exactly in the resulting potential, 
is a consistent procedure. In general, it is difficult to compare the order of accuracy 
of different approaches and further work is required in this direction. 

Our analysis uses a metric which covers not only the Schwarzschild spacetime 
but also several generalizations including Rcissner-Nordstrom, De Sitter and 
Schwarzschild-De Sitter. While the naive integration performed in Hll|l uses only 
the form of the metric near the horizon and will go through for all the cases, there 
are other subtleties which must be borne in mind while generalizing this result. In 
the case of De Sitter spacetime, for example, the lack of asymptotic flatness will be 
an issue in defining the QNM in our approach. In SdS, the situation is much worse 
since there is no definition of global temperature for the spacetime with more than 
one horizon. Naive application of (|llfl will now pick up contributions from the two 
horizons additively if we suitably restrict the range of integration in equation Hl()(l . It 
is not obvious that multi horizon spacetime can be handled in such a manner and this 
issue is under investigation [14| . 

Finally, there are (at least!) two questions of conceptual nature which we have 
not answered. Physical quantities with a quantized spectrum are of special interest 
when they have constant spacing. In the case of horizon area, for example, one can 
attempt to relate this to the intrinsic limitations in measuring length scales smaller 
than Planck Length ^2]. But the uniform spacing of QNM frequencies is a purely 
classical result and hence is harder to understand physically. While relating it to the 
poles of the Gamma function could be a step forward, we have a long way to go in 
understanding it. Second, the procedure does not provide an insight into the real part 
of k n and it does not seem easy to relate this result to previous analytic derivations 
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of Re k n . The results in ref. [S], for example, depends crucially on the form of the 
metric near r — 0, while we have only used the form of metric in the range r > a. 
I thank T. Roy Choudhury for comments on the draft of the paper. 
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